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Abstract: Expansions that furnish increasingly good approximations to real 
numbers are usually related to dynamical systems. Although comparing dy- 
namical systems seems difficult in general, Lochs was able in 1964 to relate 
the relative speed of approximation of decimal and regular continued fraction 
expansions (almost everywhere) to the quotient of the entropies of their dynam- 
ical systems. He used detailed knowledge of the continued fraction operator. 
In 2001, a generalization of Lochs' result was given by Dajani and Fieldsteel in 
, describing the rate at which the digits of one number-theoretic expansion 
determine those of another. Their proofs are based on covering arguments and 
not on the dynamics of specific maps. In this paper we give a dynamical proof 
for certain classes of transformations, and we describe explicitly the distribu- 
tion of the number of digits determined when comparing two expansions in 
integer bases. Finally, using this generalization of Lochs' result, we estimate 
the unknown entropy of certain number theoretic expansions by comparing 
the speed of convergence with that of an expansion with known entropy. 



1. Introduction 

The goal of this paper is to compare the number of digits determined in one ex- 
pansion of a real number as a function of the number of digits given in some other 
expansion. For the moment the regular continued fraction (RCF) expansion and the 
decimal expansion will serve as examples, the digits then being partial quotients and 
decimal digits. It may seem difficult at first sight to compare the decimal expan- 
sion with the continued fraction expansion, since the dynamics of these expansions 
are very different. However, in the early nincteen-sixtics G. Lochs [Tll | obtained 
a surprising and beautiful result that will serve as a prototype for the results we 
are after. In it, the number m(n) of partial quotients (or continued fraction digits) 
determined by the first n decimal digits is compared with n. Thus knowledge of 
the first n decimal digits of the real number x determines completely the value 
of exactly m(n) partial quotients; to know unambiguously more partial quotients 
requires more decimal digits. Of course m(n) — m(x, n) will depend on x. But 
asymptotically the following holds (in this paper 'almost all' statements are with 
respect to the Lebesgue measure A). 

Theorem 1.1 (Lochs). For almost all x: 
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Roughly speaking, this theorem tells us that usually around 97 partial quotients 
are determined by 100 decimal digits. By way of example Lochs calculated (on an 
early main- frame computer (J^) that the first 1000 decimals of ir determine 968 
partial quotients. 

Examining the function m(n) more closely, one can formulate Lochs' result as 
follows. Let B n {x) be the decimal cylinder of order n containing x, and C m (x) the 
continued fraction cylinder of order m containing x, then m(n) is the largest index 
for which B n (x) C C m (x). Denoting by S the decimal map S(x) = lOx mod 1, 
and by by T the continued fraction map T(x) = - — [-\ , then Lochs' theorem says 
that 

m(x) h(S) 
hm = — — (a.e.). 

n— >oo n 11(1 ) 

In 2001, the second author and Fieldsteel generalized in Q Lochs' result to any two 
sequences of interval partitions satisfying the result of Shannon-McMillan-Breiman 
theorem, (see p. 129), that can be formulated as follows. 

Theorem 1.2 (Shannon McMillan Breiman). Let T be an ergodic measure 
preserving transformation on a probability space [X, B, fx) and let P be a finite or 
countably infinite generating partition for T for which (P) < oo, where (P) 
denotes the entropy of the partition P. Then for /i-almost every x: 

hm - l0g ^ (P " (a;)) = K (T) . 
n— >oc n 

Here h^ (T) denotes the entropy ofT and P n (x) denotes the element of the partition 
Y™^ 1 T~ l P containing x. 

The proof in |7j is based on general measure-theoretic covering arguments, and 
not on the dynamics of specific maps. However, with the technique in [7| it is quite 
difficult to have a grip on the distribution of m(n). For that, one needs the dynamics 
of the underlying transformations as reflected in the way the partitions are refined 
under iterations. 

In this paper we give a dynamical proof for certain classes of transformations, 
and we describe explicitly the distribution of m(n) when comparing two expansions 
in integer bases. We end the paper with some numerical experiments estimating 
unknown entropies of certain number theoretical expansions by comparing them 
with expansions of known entropy. 



2. Fibred systems 

In order to generalize Lochs' result we start with some notations and definitions that 
place maps like the continued fraction transformation in a more general framework; 
see also |l5|, from which the following is taken with a slight modification. 

Definitions 2.1. Let B be an interval in K, and let T : B — > B be a surjective 
map. We call the pair (_B, T) a fibred system if the following conditions are satisfied: 

(a) there is a finite or countably infinite set D (called the digit set), and 

(b) there is a map k : B — > D such that the sets 

B(i) = k^ii} = {x e B : k(x) = i] 



form a partition of B; we assume that the sets B{i) are intervals, and that 
V = {B(i) : i G D} is a generating partition. Moreover, we require that 
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(c) the restriction of T to any B(i) is an injective continuous map. 

The cylinder of rank n determined by the digits k\ , k 2 , . . . , k n € D is the set 

B(kt, hi,..., k n ) = B(h) n T- x B(k 2 ) n . . . n T- n+1 B(k n ); 

by definition B is the cylinder of rank 0. Two cylinder sets B n ,B* of rank n are 
called adjacent if they are contained in the same cylinder set £? n _i of rank n — 1 
and their closures have non-empty intersection: 

For any fibred system (B, T) the map T can be viewed as a shift map in the following 
way. Consider the following correspondence 

* : (k 1} k 2 , . . . ,k n , . . .), 

where fcj = j if and only if T % ~ x x e B(j), then 

T:*(a;) !->(*&,*&,...,&„,...). 

We write B n (x) for the cylinder of rank n which contains x, i.e., 

B n (x) = B(ki, k 2 , ■ ■ ■ ,k n ) $>(x) = (k 1 ,k 2 ,...,k n ,...). 

A sequence (ki, . . . , k n ) e D n is called admissible if there exists x e B such that 
= (fci,..., /:„,... )■ 

It will be clear that ^ associates to any x G B the 'expansion' (k\, k 2 , ■ ■ ■), and 
that the expansion is determined by T and the digit map k. In order to arrive at 
the desired type of result for the fibred system (B,T), we need to impose some 
restrictions on T. 

Definition 2.2. A fibred system (B,T) is called a number theoretic fibred system 
and the map T a number theoretic fibred map if T satisfies the additional conditions: 

(d) T has an invariant probability measure \xt which is equivalent to A, i.e., there 
exist ci and c 2 , with < c± < c 2 < oo, such that 

ciX(E) < fi T {E) < c 2 \{E), 

for every Borel set E C B, and 

(e) T is ergodic with respect to Ht (and hence ergodic with respect to the 
Lebesgue measure A). 

Examples 2.3. The following examples are used throughout the text. 

(i) The first example is the regular continued fraction. Here B = [0, 1) C R, the 
map T is given by T(x) = ± - k(x), and k(x) = [jj. Hence B(i) = (j^, I], 
for all i > 1. It is a standard result that the map is ergodic with the Gauss- 
measure as invariant measure. The entropy of T equals h(T) — tt 2 / (6 log 2) = 
2.373- ■-. 

(ii) The second standard example is the g-adic expansion, for an integer g > 2. 
Here the map is T g {x) — g ■ x — k(x), and k(x) = \_g ■ x\. Hence T g (x) is the 
representative of g ■ x mod I in B = [0, 1), and B(i) = [^, ^-), for < i < g. 
Again this is a number theoretic fibred map, with invariant measure A. The 
entropy of T g equals h(T g ) = log g. 
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(iii) A third example is given by the Liiroth expansion. Now B(i) = [jrrp j) for 
i > 1. If T is taken as an increasing linear map onto [0, 1) on each B(i), 
the Liiroth series expansion is obtained; if T is taken to be decreasing linear 
onto [0, 1) it generates the alternating Liiroth series expansion of a number in 
[0, 1). See Q and H3, where it was shown that the invariant measure is the 
Lebesgue measure and that the entropy equals 

, . ^ Wfcffc + 1)) 



(iv) As a fourth example we take Bolyai's expansion; again B — [0, 1), and now 
T : [0,1) -> [0,1) is defined by 



where 



T(x) = (x+ l) 2 - l-£i(x), 

r 0, if a: e B(0) = [0, \/2- 1); 
ei(a:)=< 1, if xe = [y/2- 1,V3- 1); 
[ 2, if i£ S(2) = [V3- 1,1). 

Setting e„ = e„(x) = £i(T n_1 (:r)), n > 1, one has 



£i + y £2 + V£3 H • 

It is shown in 01 that T is a number theoretic fibred map. Neither the 
invariant measure [i nor the exact value of the entropy h(T) is known; compare 
Experiment l|6.3|l and the reference to below. 

(v) As a final example, let j3 > 1 satisfy 1 = (3 + 1, i.e., j3 is the golden mean. 
We describe the (3-continued fraction expansion of x. It is an interesting com- 
bination of regular continued fractions and the /3-adic expansion. A (3-integer 
is a real number of the form 

a n /3 n + an-xP 11 - 1 + ■ ■ ■ + a 

where a n , a n _i, . . . , ao is a finite sequence of 0's and l's without consecutive 
l's, i.e., ai ■ a,_i ^ 1 for i = 1, 2, . . . , n. In fact, /3-integers can be defined for 
any (3 > 1; see e.g. @ or 

For & S (0,1), let [^J^ denote the largest /3-integer < i, and consider the 
transformation T : [0, 1) -> [0, 1), defined by T(x) = | - [ij for i^0, and 
T(0) = 0. Iteration of T generates continued fraction expansions of the form 



b 3 + 



where the bi are /3-integers. This expansion is called the /3-continued fraction 
of x. 

Clearly, this continued fraction expansion is not an /-expansion (cf. |l5lp in 
the classical sense, but it has many of the essential properties of /-expansions, 
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and thus Renyi's approach in [14( suggests that T has a finite invariant mea- 
sure equivalent to the Lebesgue measure, with a density bounded away from 
zero and infinity. Very little seems to be known about such continued frac- 
tions. Recently, Bernat showed in Q that if x € Q(/3), then x has a finite 
/3-expansion. 

We have the following lemma. 

Lemma 2.4. Let T be a number theoretic fibred map on B, then for almost all x: 

iim = 1. 

n-oo log fl T (Bn (Zj) 

Proof. Since 

Cl \(B n {x)) < i*r{B n {x)) < c 2 X(B n (x)), 

it follows that 

-logc 2 +logfj, T {B n (x)) log\(B n (x)) - logci + log/j T (B n (x)) 



log \x T {B n [x)) \og(i T (B n (x)) log fi T (B n (x)) 

Taking limits the desired result follows, since V is a generating partition and there- 
fore lim, woo fjrr(B n (x)) = almost surely. □ 

The following formalizes a property that will turn out to be useful in comparing 
number theoretic maps. 

Definition 2.5. Let T be a number theoretic fibred map, and let / C B be an in- 
terval. Let m = m(I) > be the largest integer for which we can find an admissible 
sequence of digits oi, a-i, ... , a m such that / C B(a±, 02, ... , a m ). We say that T is 
r-regular for r E N, if for some constant L > 1 the following hold: 

(i) for every pair B, B* of adjacent cylinders of rank n: 

1 < A( jP < L 
L ~ A(5*) " ' 

(ii) for almost every x € I there exists a positive integer j < r such that either 
B m+j (x) d I or B* n+] C / for an adjacent cylinder B* m+0 of B m+J (x). 

Thus r-regularity expresses, loosely speaking, that if the expansion T agrees to 
the first to digits in both endpoints of a given interval /, then for a.e x £ I there 
exists a cylinder of rank at most 777 + r with the property that it or an adjacent 
cylinder contains x, and is contained entirely in /; moreover, the size of two adjacent 
cylinders differs by not more than the constant factor L. 

Examples 2.6. We give three examples. 

(i) Let / C [0, 1) be a subinterval of positive length, and let m = m{I) be such 
that B m is the smallest RCF-cylinder containing /, i.e., there exists a vector 
{a\, . . . , a m ) e N m such that B m = B m (a\, . . . , a m ) D I and A(J n ([0, 1) \ 
i? m _l_i(ai, . . . , a m , a))) > for all a € N. One can easily check that r = 3 in 
case of the RCF. 

Furthermore, it is well-known, see e.g. (4.10) in p. 43, that 
X(B m (ai, . . . , a m )) 



Qm(Qm Qm—x) 
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which yields, together with the well-known recurrence relations for the partial 
fraction denominator sequence (Q m )m>o (see 0], (4.2), p. 41): 



that 



-l 0; Qq :— 1; Q m — a m Q m -i + Q m -2, m > 1, 



X(B m (ai, . . . , a m )) < 3X(B m (ai, ...,a m + 1)), 



for all (<xi, . . . , a m ) £ N m , i.e., the continued fraction map T is 3-regular, with 
L = 3. 

(ii) The g-adic expansion is not r-regular for any r. Although adjacent rank n 
cylinders are of the same size, so L = 1 can be taken in I j2.5|l (i). property 
(|2.5(l fii'l does not hold. This can be seen by taking a; in a very small interval 
/ around 1/g. 

(iii) The alternating Liiroth map is 3-regular with L = 2. 
3. A comparison result 

We would like to compare two expansions. In this section we will therefore assume 
that S and T are both number theoretic fibred maps on B. We denote the cylinders 
of rank n of S by A n , and those of T by B n . We let m(x, n) for x £ B and n > 1 be 
the number of T-digits determined by n digits with respect to S, so m(x, n) is the 
largest positive integer m such that A n (x) C B m (x). Yet another way of putting 
this, is that both endpoints of A n (x) agree to exactly the first m digits with respect 
to T. 

We have the following general theorem. 
Theorem 3.1. Suppose that T is r-regular. Then for almost all x € B 

m(x,n) h(S) 



lim 



h(TY 



Proof. For x £ B let the first n S'-digits be given; then by definition of m = m(x, n) 
one has that A n (x) C B m (x). Since T is r-regular, there exists 1 < j < r such that 

B m+j C A n (x) C B rn (x), 

where B m+ j is either B m+ j{x), or adjacent to it of the same rank. By regularity 
again, for some L > 1, 

-\{B m+r (x)) < -X(B m+j {x)) < X(B m+j ), 

so 

-(-lo gj L + logA(B m+r (x))) < -logA(A„(x)) < - \ogX(B m (x)), 
n n n 

and the result follows from Lemma (|2.4|) (applied to both T and S), and the The- 
orem of Shannon-McMillan-Breiman 11. 2|) . □ 

Example 3.2. Let m^ CF (x, n) be the number of partial quotients of x determined 
by the first n digits of x in its <?-adic expansion. Then for almost all x: 



lim 



m£ CF (x,n) _ 6 log 2 log g 
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This generalizes Lochs' theorem to arbitrary g-adic expansions. 

Example 3.3. Let m g IjE (x,n) be the number of alternating Liiroth digits of x 
determined by the first n digits of x in its g-adic expansion. Then for almost all x: 



lim 



m£ LE (x,n) log g 



fc=l 



3gfc(fc+l) 
fc(fc+l) 



Note that the conditions in Theorem l3. H are not symmetric in S and T. Since the 
regularity condition is not satisfied by /i-adic expansions, Theorem 13. II is of no use 
in comparing g-adic and /i-adic expansions, nor in proving that the first n regular 
partial quotients determine usually 7r 2 /(61og21ogg) digits in base g. 

4. Comparing radix expansions 

In this section we compare g-adic and /i-adic expansions, for integers g, h > 2, by 
explicitly studying the distribution of m(n). 

Given a positive integer n, there exists a unique positive integer t — £(n), such 
that 

h~ {e+1) < g- n < h- e . (4.1) 

Thus the measure X(A n ) of a g-cylinder of rank n is comparable to that of an 
/i-cylinder of rank £(n). 
It follows that 

£{n) log 5 • t i(n) + l 
lim sup < < lim mt , 

n^oo n log h n— >oo n 

i.e., 

hm = -, (4.2) 

n log/i 



11 — >OQ 



which is the ratio of the entropies of the maps T g and Th introduced in l|2.3|l (ii) . 
Of course, one expects the following to hold: 

Um 4^M = l£i| (a . e .). (4.3) 
n^oo n log h 

which is the analog of Lochs' result for the maps T g and T^; here m g h \x,ri) is 
defined as before: it is the largest positive integer m such that A n (x) is contained 
in the /i-adic cylinder B rn (x). One should realize that, in general, £(n) has no 
obvious relation to m(n); in fact equation l|4.2Jl is merely a statement about the 
relative speed with which g-adic and /i-adic cylinders shrink. 

Let x G [0, 1) be a generic number for S = T g for which we are given the first n 
digits ti,t%, . . . ,t n of its g-adic expansion. These digits define the g-adic cylinder 

A n (x) = A n (t 1 ,t 2 , ■ ■ ■ ,i„). 

Let m(n) = m g (n) and l(n) be as defined above; note that 

m(n) < £(n), for all n > 1. (4.4) 

The sequence (m(n))^° =1 is non-decreasing, but may remain constant some time; 
this means that it 'hangs' for a while, so m(n + t) = ■ ■ ■ = m(n + 1) = m(n), after 
which it 'jumps' to a larger value, so m(n + 1 + 1) > m(n + t). Let (nk)k>i be the 
subsequence of n for which m(n) 'jumps', that is, for which m(n^) > m(rifc — 1). 
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Lemma 4.5. For almost all x, 

m { g l \n k ) log g 
km = - — -. 

k~ »oo n k log h 

Proof. By definition, £> m („ fc )(x) D A rik (x). The cylinder B m ^ nk ^(x) consists of h 
cylinders of rank m{rik) + 1, and A nk (x) intersects (at least) two of these, since 
otherwise some -B m (n fc )+i D A nk (x), contradicting maximality of mink), and thus 
it contains an endpoint e of some B m / nk \ + i lying in the interior of B m ^ nk ^(x). 

On the other hand, by definition of n k , we know that m{n k ) > m(n k — 1), 
so B m ( nk }(x) 7$ A rik _i(x), and therefore A nk ^i(x) contains an endpoint / of 
B m (n k ){x) as well as e. Now e and / are at least \{B m i nh \ + x) = ft,~( m ( ,lfc ) +1 ) apart. 
Therefore 

h -(m(n k )+i) < X (A nk ^(x)) = g ■ X(A nk (x)) = g ■ g- n * < g ■ h~ l ^\ 
by (|4.1|l . Hence 

^-(m(n*.)+l) < _ fo-l(n k ) 

which, in combination with (|4.4|) implies 

*(nk) - 1 - < m(n*) < i(n fc ), k > 1. (4.6) 

But from l|4.6|) the Lemma follows immediately. □ 

Next we would like to show that l|4. 3JI follows from 1)4. ti [I : this is easy in case 
h = 2. Compare Q. 

Corollary 4.7. W*£/i notations as above, for any g £ N, g > 2 

Hm mgW0 = l2M (a . e .). 
n—>oo n log 2 

Proof. If ft, = 2 the mid-point £ of -B m ( n ) lies somewhere in A„. Now m(n + 1) = 
m(n) if this mid-point £ is located in A n +i(ti, ti, . . . ,t n , £ n +i), that is, if it is located 
in the same g-adic cylinder of order n + 1 as i. Notice that this happens with 
probability |, and that the randomness here is determined by the <?-adic digit t n+ i. 
To be more precise, let H denote the event that we will 'hang' at time n, i.e., the 
event that m(n) — m(n + 1), and let D be a random variable with realizations 
t € {0,1,..., 3- 1}, defined by £ <E A n+1 (tx, . . . ,t n , D). Now 

9-1 

P(H)=Y,?{B\D = i)'P(D = i), 

i=0 

and from P(H\D = i) = l/g for < i < g - 1 it then follows that P(H) = l/g; due 
to the discrete uniform distribution of the digit t n+ i of x we do not have to know 
the probabilities P(D = i). 

Thus we see that for each k > 1 (and with no = 0) the random variable 

Vk = n k - Uk-i 

is geometrically distributed with parameter p = l/g. Furthermore, the v^s are 
independent, and, since n k > k 

n k+ i n k + n k+ i - n k 1 1 

1 < = = H Wfc+l < 1 + 7-Wfc+l- 

n k n k n k k 
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Because V\ , V2 , ■ ■ ■ are independent and identically distributed with finite expecta- 
tion, it follows from the Lemma of Borel-Cantelli that lim^oo -j^ 1 - = (a.e.), and 
therefore 

v Uk+1 1 1 \ 
lim = 1 (a.e. J. 

k — >oo n k 

Given any n > 1, there exists k = k(n) > such that n k < n < n k+ \. Since 
fn{n k ) < m(n) < m{n k+ i) one has 

m(n k ( n )) m(n) m(n k{n)+1 ) 
nk(n)+i ~ n ~ 7ife(„) 

But then the result follows from (|4. f>|> with ft = 2. □ 

If ft. > 3 the situation is more complicated; there might be more than one mid- 
point £i 'hitting' A n . In case only one ^ 'hits' A n (as is always the case when ft, = 2), 
we speak of a type 1 situation. In case there is more than one 'hit' we speak of a 
type 2 situation. Now change the sequence of jump-times (n k ) k >i by adding those 
n for which we are in a type 2 situation, and remove all n k s for which we are in 
a type 1 situation. Denote this sequence by (n k ) k >i. If this sequence is finite, we 
were originally in a type 1 situation for n sufficiently large. In case (n k ) k >i is an 
infinite sequence, notice that for every n for which there exists a k such that n = n* k 
one has 1-regularity, i.e., 

B m (n' k ) + 1 C A n *(x) C B m ( n *)(x) , 

where £? m ( n *)+i is either B m / n *\+i(x) or an adjacent interval. Following the proof 
of Theorem fo.ll one has 

mint) logo , , 
hm — = (a.e.). 
fc— »oo n k log h 

Let {hk)k>i be the sequence we get by merging (n k )k>\ and {n* k )k>i- Notice that 

mini.) logo , , . . 

lim — )r^- = — (a.e.). (4.8) 

k^oo n k log ft 

If {hk)k>i is N \ (hk)k>i, we are left to show that 

m(h k ) \ogg 
hm — - = (a.e.). 

k^oo n k log ft 

Let n = hj for some j > 1. Then there exist unique k = k{j) and ft = ft(j') such 
that n k+ i = n/i+i, n fe < n/i and 

n h <n < n k+1 . 

Since n k+ i — hh = hh+i — hh = Vh+i is geometrically distributed with parameter 
1/g, and since v\, V2, ■ . ■ are independent with finite expectation, it follows as in the 
case ft = 2, that 

— — = 1 + —v h +i < 1 + t u 'i+i — * 1 as j — » 00 (a.e.). 
n/j n ft ft 

But then (|4.3|l follows from 14. 8|) and from 

m(h h ) m(n) m(h h+1 ) 
hh+i ~ n ~ hh 

We have proved the following theorem; again, see also 0. 
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Theorem 4.9. Let g,h e N>2- TTien /or almost all x 

rn g h \x,n) logg 
hm = -. 

n—too n log h 

5. Remarks on generalizations 

At first sight one might think that the proof of Theorem l4.9l can easily be extended 
to more general number theoretic fibered maps, by understanding the distribution of 
m(n). However, closer examination of the proof of Theorem l4.9l shows that there are 
two points that make generalizations of the proof hard, if not impossible. The first is 
the observation, that the 'hanging time' Vk (in a type 1 situation) is geometrically 
distributed, and that the v^s are independently identically distributed. Recall that 
this observation follows from the fact that the digits given by the g-adic map S g 
are independently identically distributed and have a discrete uniform distribution. 
As soon as this last property no longer holds (e.g., if S = T 7 , the expansion with 
respect to some non-integer 7 > 1), we need to know P(D = i) for each i S 
{0, 1, . . . , [7] — 1}, and this might be difficult. 

Even if we assume S to be the g-adic map, another problem arises when T is 
not the /i-adic map, but, for example, T — T 7 . In that case all the ingredients of 
the proof of Theorem 14 . 91 seem to work with the exception that at the 'jump-times' 
(n-k)k we might not be able to show that 

t{n k ) - C < m(n k ) < £(n k ), (5.1) 

with C < 1 some fixed constant. Notice that from (|5.1() it would follow that 

y rn(n k ) _ logg 
fe^oo n k log 7 

For one class of 7 € (1,2) one can show that (|5.1(l still holds. These 7's are the 
so-called 'pseudo-golden mean' numbers; 7 > 1 is a 'pseudo-golden mean' number 
if 7 is the positive root of X k — X k_1 — ■ ■ • — X - 1 = 0, for some k E N, k > 2 
(in case k = 2 one has that 7 = (3 = |(1 + \/5), which is the golden mean). These 
'pseudo-golden mean' numbers 7 are all Pisot numbers, and satisfy 

, 1 1 1 
1 = - + — + ■■■ + — • 

Theorem 5.2. With notations as before, let S = T g and T — T 7 , with g G N>2, 
and 7 > 1 a 'pseudo golden number'. Then 

r4 7) (n) logg 
lim = , (a.e.). 

n^oo n log 7 

This result is immediate from pj- Here we present a different proof, based on 
the distribution of the number of correct digits, m(n). 

Proof of Theorem \5. l <\ In case k = 2, a digit 1 is always followed by a zero, and 
only the 7-cylinders Bm' corresponding to a sequence of 7-digits ending with are 
refined. Thus from the definition of m(n) it follows that the last digit of B^}jx) 
is always (if it were 1, the choice was wrong). Note that in this case 



A(^(x))=7- m( " ) . 
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Let £(n) be defined as before, and notice that at a 'jump-time' n = rih one has 

7 -Mn)+2) <g.g-n <g. 7 -*( n ), (5.3) 

from which 

£( n ) _2 - < m(n). 
log 7 

Since m{n) < £{n) we see that (|5.1|) follows in case 7 equals the golden mean. 

In case k > 3 the situation is slightly more complicated; let us consider here 
k = 3. Now any sequence of 7-digits ending with two consecutive l's must be 
followed by a zero, and therefore — by the definition of m(n) — the last digit of 
or is 1 which is preceded by 0. One can easily convince oneself 

then that 

XCrW Ml - / 7 ~ m(n) if d Mn)( x ) = °> 

1 ™(«) 7 -(™(«)+l) + 7 "(m(n)+2) jf dm{n) ( x ) = l. 

Here d m („)(x) is the m(n) th 7-digit of x. 

Now let n — be a 'jump-time'. If d TO ( n )(x) = 0, then B^.Ax) consists of two 
7-cylinders, one of length T -Mn)+l) an d one of length 7 ~M*)+ 2 ) + 7 -(™(™)+3)_ 
One of these two cylinders is contained in A n -\(x) (due to the fact that n is a 
'jump-time'), and (|5.3|l is therefore satisfied. Of course one has m(n) < £(n) in this 
case. In case d m („) (x) = lwe see that B 1 ^,-. (x) consists of two 7-cylinders, one of 
length 7 -M«)+i) an d one of length 7 -(™(™)+2). Now 

A(B (7 ( (x)) = 7 -( m ( rl ) +1 ) + 7 -(™(™)+2) < 7 -m(n) j 

and by definition of £{n) one has m(n) < £(n). Again one of these two sub-cylinders 
of B%l n) (x) is contained in A„_i(i), and one has 

7 -(m(n)+2) < g ■ g~ n < g ■ 7 -*(- n l 



We see that (|5.3|) is again satisfied. For k > 4 the proof is similar; one only needs 
to consider more cases. □ 



6. Estimation of (unknown) entropies 

In this section we report on some numerical experiments. Our experiments were 
carried out using the computer algebra system Magma, see . The general set-up 
of the experiments was as follows. We choose n random digits for the expansion of 
a number x in [0, 1) with respect to a number theoretic fibred map S, and compute 
m'g(x,n), the number of digits of x with respect to the number theoretic fibred 
map T, determined completely by the first n digits of x with respect to S. This is 
done by comparing the T-expansions of both endpoints of the S'-cylinder A n (x). 

Experiment 6.1. The first experiment was designed to test the set-up. It aimed 
to verify Lochs' result: we let S = T\q, the decimal expansion map, and T = T RCF 
for the regular continued fraction map. For N = 1000 random real numbers of 
n = 1000 decimal digits, we computed the number of partial quotients determined 
by the endpoints of the decimal cylinder A n {x). 

We ran this experiment twice; the averages (over 1000 runs) were 970.534 and 
969.178 (with standard deviations around 24). Compare these averages to the value 
970.270 • • • predicted by Theorem ITU 



Entropy quotients and correct digits 
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Experiment 6.2. In this experiment we compared the relative speed of approxi- 
mation of g-adic and h-adic expansions: 



h -> 




2 




7 




10 


9l 


pred 


observed 


pred 


observed 


pred 


observed 


2 






0.356 


0.355 0.355 


0.301 


0.300 0.300 


7 


2.807 


2.805 2.805 






0.845 


0.844 0.844 


10 


3.322 


3.320 3.320 


1.183 


1.182 1.182 







In fact we compared binary, 7-adic and decimal expansions of N — 1000 random 
real numbers with each of the other two expansions, again by computing to how 
many /i-digits both endpoints of the g-adic cylinder Aiqqq(x) agreed. The two tables 
list the values found in each of two rounds of this experiment and compares the 
result with the value found by Theorem 14.91 Note that by (|4.4I) always m(n) < 
£(n) and thus the observed ratio should approximate log g/ log h from below. By 
reversing the roles of g and h, it is thus possible in this case to approximate the 
entropy quotient by sandwiching! 

Experiment 6.3. By choosing random real numbers of 1000 decimal digits, we 
attempted to estimate the unknown entropy of the Bolyai map (cf. (|2.3|l fiv)'). We 
found that on average 1000 digits determined 2178.3 Bolyai digits in the first run 
of N — 250 random reals, and 2178.0 in the second; the standard deviations were 
around 13. These values would correspond to an entropy of around 1.0570 or 1.0572. 

Choosing N random Bolyai expansions of length 1000 (by finding the Bolyai 
expansion to 1000 digits for random reals), we determined the number of decimal 
digits in one experiment, and the number of regular continued fraction partial 
quotients in another experiment. In the first run of N = 1000, we found that on 
average 458 decimal digits were determined, in the second run 457.8 on average 
(with standard deviations around 4) . These values would correspond to an entropy 
of around 1.054G for the first, and 1.0541 for the second run. according to Theorem 
5 in 0. In runs with TV = 1000 of the second experiment we found that 444.36 and 
444.37 partial quotients were determined on average, with standard deviations of 
around 17. This suggests an entropy of 1.0545 or 1.0546 according to Theorem 13. II 

These entropy estimates for the Bolyai-map were reported in a preliminary ver- 
sion of this paper |5(. Using a different method, based on the fixed-points of the 
Bolyai-map, Jenkinson and Pollicott were able to find in [£j sharper estimates for 
the entropy. To be precise, they found 1.056313074, and showed that this is within 
1.3 x 10~ 7 of the true entropy. 

Experiment 6.4. In our final experiment we attempted to estimate the entropy 
of the /3-continued fraction map, see (|2.3li fv~) . 

As before, we chose real numbers with 1000 random decimal digits and deter- 
mined the number of /3-integer partial quotients determined by them. In the first 
run, the 1000 decimals determined 877.922 partial quotients on average, in the sec- 
ond run 878.125 partial quotients were determined. These values correspond to an 
entropy of 2.021 and 2.022. 
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